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ABSTRACT 

Spectroscopy of the low mass X-ray binary SerX-1 using the Gran Telescopio Canarias 
have revealed a ~2 hr periodic variability that is present in the three strongest emission 
lines. We tentatively interpret this variability as due to orbital motion, making it the 
first indication of the orbital period of Ser X-1. Together with the fact that the emission 
lines are remarkably narrow, but still resolved, we show that a main sequence K-dwarf 
together with a canonical IAMq neutron star gives a good description of the system. 
In this scenario the most likely place for the emission lines to arise is the accretion 
disk, instead of a localized region in the binary (such as the irradiated surface or 
the stream-impact point), and their narrowness is due instead to the low inclination 
(s^lO°) of Ser X-1. 
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1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are interacting binaries 
where a low-mass donor transfers matter onto the neutron 
star or black hole via Roche-lobe overflow. This process 
forms an accretion disk and gives rise to the observed X- 
rays, making them among the brightest X-ray sources in 
the sky (see e.g. the collected reviews in Lewin & van der 
Klis 2006). Although optical counterparts for many of the 
bright LMXBs have been known for a long time, kinematic 
studies of these systems are often not straightforward. This 
is mainly due to reprocessing of the X-rays in the outer 
accretion disk that completely dominates the optical light. 
This situation has changed only with the advent of sensitive 
spectrographs on very large telescopes. 

Using phase-resolved medium resolution spectroscopy 
of ScoX-1, Steeghs & Casares (2002) detected the pres- 
ence of narrow high excitation lines that were strongest in 
the Bowen region (a blend of NIII A4634/4640 and CUT 
A4647/4650). These narrow lines originate from the irradi- 
ated surface of the donor star and allowed for the first time 
a radial velocity study of both binary components in Sco X- 
1 which led to constrain its mass function. Phase-resolved 
spectroscopic studies of about a dozen other bright LMXBs 
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have revealed the presence of these narrow components in 
the Bowen region in most LMXBs and provided in many 
cases the first constraints on the mass of their compact ob- 
jects (e.g. Cornelisse et al. 2008 for an overview). 

One such bright LMXB is SerpensX-1 (Ser X-1). It was 
discovered in the 1960s (Friedman et al. 1967), but its opti- 
cal counterpart (MM Ser) was only identified about 10 years 
later (Davidsen 1975). Further observations taken in good 
seeing conditions showed that the optical counterpart con- 
sisted of a blend of two stars that are separated by ~2.1 
arcsec (Thorstensen et al. 1980). Subsequently, one of these 
was itself revealed to be a blend of two stars that are only 
separated by 1 arcsec (Wachter 1997). Both spectroscopy 
(Hynes et al. 2004) and radio observations (Migliari et al. 
2004) confirmed that the bluer of the stars is the true optical 
counterpart. Due to the well-established presence of ther- 
monuclear X-ray bursts from Ser X-1, the compact object is 
known to be a neutron star (Li & Lewin 1976). Furthermore, 
Ser X-1 was one of the first neutron star LMXBs in which 
a broad relativistic iron line was detected, and modeling of 
this line inferred an inclination of <25° (Bhattacharyya & 
Strohmayer 2007; Ng et al. 2010). However, despite over 40 
years of study little else is known about Ser X-1. 

Here we present the results of a spectroscopic study of 
Ser X-1 using the Gran Telescopio Canarias (GTC) at the 
Observatorio del Roque de los Muchachos on the Canary 
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Figure 1. Average normalized OSIRIS/GTC spectrum of SerX-1. Indicated are the most prominent identified emission lines, sky lines 
(I/S) and diffuse interstellar bands (DIB) while the unidentified features are labeled with "?" . We note that several of these unidentified 
features were previously observed by Hynes et al. (2004). 



Island of La Palma (Spain). We detect a ~2 hr periodic 
variation that we interpret as the orbital period. Although 
the emission lines are very narrow and sharp (as was already 
reported by Hynes et al. 2004), we present further evidence 
that this is due to the low inclination of the system instead 
of the lines being formed on the irradiated surface of the 
donor star or the stream- impact point. 



2 OBSERVATIONS AND DATA REDUCTION 

On 6 July 2011 we observed SerX-1 for 4 hrs continuously 
(UT 21:45-02:23) using OSIRIS on GTC. For each exposure 
we used the volume-phased holographic grism R1200V with 
a slit width of 1.0 arcsec, resulting in a total of 24 spectra 
with an exposure time of 600 seconds each. The seeing dur- 
ing the observations varied between 0.7 and 1.0 arc-seconds. 
Since SerX-1 is ~1 arcsec from a fainter field star, we aligned 
the slit in such a way that we minimized the contribution of 
this interloper. However, this does mean that no comparison 
star is included in the slit and we could therefore not correct 
for slit losses. Arc lamp exposures were taken for wavelength 
calibration during the daytime. 

For the data reduction (i.e. debias, flatfielding etc) we 
used the PAMELA software that allows for optimal extrac- 
tion of the spectra (Home 1986). We determined the pixel- 
to-wavelength scale using a 4th order polynomial fit to 20 
lines resulting in a dispersion of 0.79 A pixel" ^, an rms scat- 
ter of <0.05 A, and a wavelength coverage of AA4434-6060. 
The final resolution around the Hell A4686 emission line 
is 107.6 km s~^ (FWHM). Finally we also corrected for any 
velocity shifts due to instrument fiexure (always <6 km s~^) 
by cross-correlating the sky spectra. For the corresponding 
analysis of the resulting dataset we used the MOLLY pack- 
age. 



3 DATA ANALYSIS 

3.1 Spectrum 

First we created an average normalized spectrum of SerX-1, 
and show the result in Fig.[T] The spectrum is very similar 
to that presented by Hynes et al. (2004), but there are also 
a few differences. To start with the similarities, the most 
prominent emission lines that are commonly detected in 
LMXBs are also present in our spectrum (and indicated in 
Fig-H]), and all of them are extremely narrow as shown in Ta- 
ble[T] In particular a close-up of the Bowen region, presented 
in Fig.[2l suggests that the individual NIII components (i.e. 
4634 and 4640A) are clearly separated, while there is also no 
hint of any C HI lines. Furthermore, we detect He I A4471 in 
absorption, and there are hints of features at the location of 
other He I transitions in our spectrum (which are indicated 
in Fig.[T]). However, since He I A4471 is on the edge of the 
wavelength coverage where the sensitivity starts to drop sig- 
nificantly, it is not possible to conclude anything else from 
this line. Finally, just like Hynes et al. (2004) we also de- 
tect the broad absorption trough just redward of H/3 and 
the unidentified emission features at 4604 and 4945 A. This 
confirms beyond doubt that they are not artefacts and must 
be true features of the system, and we have indicated them 
with a "?" in Fig.[Tl 

The two main differences between our spectrum and 
that of Hynes et al. (2004) are a much more complex struc- 
ture around H/3 and several new features (which we have 
also indicated with a "?"). The detailed structure around 
H/3 (ranging from 4800-4950 A) appears to form an inverted 
P Cygni profile, as if it were a single feature. However, within 
the current framework of LMXBs it is difficult to understand 
how a significant infall component above the disk of Ser X- 
1 at velocities well over 4000 km s~^ (which are needed to 
produce such a profile) could occur. Furthermore, from the 
He I A4471 absorption line we already know that at least 
He I A4922 will be present in absorption and contributing to 
this feature. We therefore think it more likely that the full 
complex around H/3 is made up of many blended emission 
and absorption lines. 

Of the new features that are present in our spectrum. 
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Figure 2. Close-up of the Bowen blend of SerX-1 with the posi- 
tion of the expected N III and C III lines indicated. Note that the 
spectrum has not been corrected for the systemic velocity of ~92 
km s~^. 



the broad emission components just redward of Hell A5411 
and the narrow one around 5004 A, are particularly interest- 
ing. They should have been clearly detected by Hynes et al. 
(2004) and to our knowledge have no counterpart in any 
other LMXB. To identify the narrow features at 4605 and 
5004 A, we used the atomic line list by van Hoof & Verner 
(1998) and the applied off-set of 92 km s"^ (see Sect. 3.2) 
from the HeH A4686 rest-frame wavelength as an estimate 
of the systemic velocity of SerX-1. The 4605 A feature cor- 
responds to a NHI transition (4603.8 A), while that at 5004 
A is close to a NH one (5002.7 A). We therefore tentatively 
conclude that these unknown transitions are due to ionised 
Nitrogen. 

We examined the broad feature at 5440A in the indi- 
vidual spectra by eye and note that it is extremely variable. 
Although present in all spectra, its strength and width shows 
large changes from spectrum to spectrum. We illustrate this 
in Fig.|3]by showing the close-up of this region in two differ- 
ent spectra. Although we have no explanation for the strong 
variability displayed, the profile in the lower spectrum in 
Fig.[3]does suggest that it consists of several unresolved lines. 
Indeed we note that a large number of N II/N HI transitions 
occur in this region that are close to the individual peaks 
displayed, but the quality and resolution of the individual 
spectra does not allow us to unambiguously conclude this. 



3.2 Radial velocities 

Following the report of Hynes et al. (2004) of a small drift of 
Hell A4686 over the period of ~an hour, we inspected the 
trail of this line by eye. We could clearly see some movement 
over time and therefore averaged two consecutive spectra 
together to improve the signal-to-noise, thereby creating 12 
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Figure 3. Close-up of the Hell A5411 emission line for two spec- 
tra of SerX-1. Although Hell is clearly detected in both spectra 
(and labeled as such), a broad and highly variable unidentified 
feature is also present (labeled with "?"). 



spectra in total. For each of these spectra we then deter- 
mined the full width at half- maximum (FWHM) , equivalent 
width (EW) and radial velocity. Our FWHM measurements 
show that the line is resolved and we therefore corrected our 
measurements for instrumental broadening using a nearby 
line in our arc spectrum. Since both the FWHM and EW do 
not change significantly during our observations we present 
their average value in Table[T] Finally, we note that the aver- 
age radial velocity is consistent with that observed by Hynes 
et al. (2004), and in Fig.|4]we show the resulting radial ve- 
locity curve. 

Interestingly, the radial velocity curve does appear to 
show periodic motion. We therefore checked other strong 
emission lines (i.e. the Bowen region and Hell A5411) to 
see if similar variability is present. The Bowen region also 
shows that there are no strong changes between the spectra 
in both strength and width of the individual components 
that make up the blend. We also note that there are clear 
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Table 1. Emission line properties in SerX-1 for both lines that 
are typically observed in LMXBs (known lines) and emission lines 
that are unique to SerX-1 (unknown lines). Included are the best 
fit period (P) and semi-amplitude {K) obtained via sine-wave 
fits for the lines where variability is observed. For the lines that 
are part of a blend we have also indicated the properties of the 
full region. All FWHM measurements are the intrinsic broadening 
after deconvolving the instrumental resolution. 



Line 


EW 


FWHM 


P 


K 




(A) 


(km s-i) 


(hrs) 


(km s-i) 


known lines 










Hell A4686 


1.20±0.07 


180±10 


1.98±0.13 


38±5 


NIII A4640 


0.81±0.06 


210±12 


2.19±0.15 


31±5 


Hell A54111 


1.01±0.04 


305±44 


1.98±0.20 


34±8 


unidentified lines 








A4604 


0.68±0.08 


418±62 






A4945 


0.57±0.06 


1061±94 






A5004 


0.30±0.04 


246±34 


1.84±0.342 


23±122 


A5290 


0.42±0.05 


823±85 






blends 










Bowen 


1.40±0.09 


717±57 






HeH A54113 


3.62±0.08 


2220±138 







^EW and FWHM are for the Hell emission line only. 

^Note that the sine-wave fit is formally not significant. 

^EW and FWHM are for the full structure between 5400-5450A 



velocity shifts in N III A4640, the strongest component of the 
Bowen region, and therefore measured the FWHM, EW and 
radial velocity for this Une. We present the results in Fig.|4] 
and Table[T] Finally we also measured the average EW and 
FWHM of the full Bowen blend, and Ust these results in 
Tablen 

As aheady pointed out in Sect. 3.1, the region around 
HeH A5411 is more complex due to the strongly variable 
component just redward (see Fig. [3]). However, HeH is itself 
clearly visible in the individual spectra and again it is mov- 
ing, and in Fig|4]we give its radial velocity curve. However, 
the FWHM and in particular EW measurements of the indi- 
vidual spectra gives such a large spread (see Fig.|3} that we 
have opted to only measure these values from the average 
of all spectra, giving the results in Table[T] Finally we also 
measured the FWHM and EW of the full complex between 
5400-5450 A, and these are also included in Table[T] 

We also searched for variability in the unidentified emis- 
sion lines in Fig. [T] (those indicated with "?"). Most lines are 
either too weak or show no variability, and we therefore only 
present their FWHM and EW in Table[T] The emission line 
at 5004 A is the only one that shows some evidence of vari- 
ability although at a low significance of ~2(t. Interestingly, 
a sine- wave fit to its radial velocity curve gives a period and 
amplitude that are consistent with the lines discussed above. 
For completeness we list the fit results in Table[Tl but will 
not include them in any further analysis. 

FigS] shows that the three emission lines all have a vari- 
ability with a period, semi-amplitude and phasing that are 
similar within the errors. To verify that the variability does 
not have an instrumental origin we also measured the radial 
velocities of a diffuse interstellar band (DIB) around 5780 
A. This is included as the bottom panel of Fig|4l and shows 



clearly no velocity shifts in the DIB feature to a limit of 
<9.4 km s"^ (99% confidence). Finally, for the HeIIA4686 
radial velocity curve a x'^-test indicates that the variability 
is highly significant, since fitting it to a constant velocity 
gives xS=5.7 (for 11 degrees of freedom). 

We therefore conclude that the variability in the 3 emis- 
sion lines is real, and fitted the radial velocity curve for each 
emission line with a sine wave (resulting in x^=0.85 for He II 
A4686). We list the resulting periods and semi-amplitudes in 
Table[T] and note that apart from the systemic velocity, 7, 
all parameters are similar within the errors. The differences 
in 7 are most likely due to the fact that both NIII A4640 
and Hell A5411 are blended with nearby lines, thereby pro- 
ducing an off-set. This is further strengthened by the fact 
that their FWHM is much broader than Hell A4686. Since 
we do not expect Hell A4686 to be blended we think that 
this line provides the best estimate of the systemic velocity, 
which we measure as 7=92±4 km s~^. 

Finally, we phase-folded both the Bowen blend and 
Hell A4686 in 8 bins to check for other emission features 
moving around on a similar period (but at a different phase 
and semi-amplitude). A visual inspection does not show any 
other moving component. Furthermore, we created averages 
by co-adding spectra spanning 0.25 orbital cycles around the 
expected minimum and maximum of both N HI A4640 and 
Hell A4686. At both extrema the lines can be very well de- 
scribed by a single Gaussian and the addition of a second 
Gaussian (either as a broader underlying feature or a second 
narrow peak) does not significantly improve the fit (accord- 
ing to an F-test). We therefore conclude that the emission 
lines can be represented by a single component. 



4 DISCUSSION 

We have presented evidence for a ~2 hrs periodic radial 
velocity variation in the emission lines from SerX-1 from 
our GTC observations. Although this variation is based only 
on approximately 2 cycles worth of data, and therefore still 
needs confirmation, we tentatively conclude that we have 
detected the orbital period of SerX-1. For the rest of the 
discussion we will assume that the variability is truly due 
to orbital motion and consider what the implications are for 
the properties of this system. 

Another important observation is the fact that most 
emission lines are narrow and are best described by a single 
component. This strongly suggests that they are located in 
a single region somewhere in the binary. Hynes et al. (2004) 
have already listed several potential locations: the outer ac- 
cretion disk, the irradiated surface of the donor star or the 
stream-impact point. Based on our proposed orbital period, 
we can now explore these different suggestions to find out 
what the consequences for the system parameters of SerX-1 
will be. 

Could the emission lines arise from the irradiated sur- 
face of the donor star? Narrow components arising from 
the irradiated donor star surface have been observed in 
many persistent LMXBs (see e.g. Cornelisse et al. 2008 
for an overview). These narrow components are in general 
strongest in the Bowen region, but there are also examples 
(e.g. EXO 0748—676) where the Hell lines show a large com- 
ponent arising from the irradiated surface (Mufioz-Darias 
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Figure 4. Radial velocity curves for the 3 strongest emission 
lines and a diffuse interstellar band (around 5780A) present in 
the spectrum of SerX-1. The emission lines are from top to bot- 
tom: Hell A4686, NIII A4640 and Hell A5411. Overplotted is the 
best fitting sine curve for each line. The radial velocity curve of 
the diffuse interstellar band on the other hand shows no obvious 
velocity movement. 



et al. 2009). Although the narrow components in many 
LMXBs have a FWHM comparable to that observed for 
SerX-1, their radial velocity semi-amplitudes are typically 
an order of magnitude larger (e.g. Casares et al. 2006, Cor- 
nelisse et al. 2007, Barnes et al. 2007). If the lines arise 
from the irradiated surface of the donor star, then the ob- 
served FWHM of the emission lines of ~180 km s~^ must 
be due to rotational broadening only (since there is no evi- 
dence that the lines consist of multiple components) . In this 
case we can easily show that for any reasonable mass ratio 
of SerX-1, the FWHM is incompatible with the observed 
radial velocity semi-amplitude of Kein=38 km s~^. 

As a first order estimate to interpret our measurements, 
we must assume a large mass ratio and that the accretion 
disk in SerX-1 is not shadowing the donor star, i.e. Kcm 
is tracing emission coming from LI. This will lead to the 
largest so-called A'-correction in order to derive the true 
radial velocity semi-amplitude of the donor star. Using the 
7^-correction polynomials by Mufioz-Darias et al. (2005) for 
g=0.8 (the largest mass ratio for which their iC-correction 



is still valid to 1%) gives a radial velocity semi-amplitude of 
the donor star of 7^2 = 111 km s~^. From this value of K2, 
together with q=0.8 we can estimate a rotational broadening 
of the emission hues of Vsim= 0.462iC2g^/^(l+g)^/^=70 km 
s~^ (Wade & Home 1988). This is much smaller than the 
180 km s~^ we observe, and suggests that g3>0.8. Even for 
the highly unlikely assumption of g~1.0 these numbers will 
still not work, and we cannot obtain a sensible solution, and 
we therefore reject the suggestion that the lines arise from 
the donor star. 

If we instead assume that the emission lines are pro- 
duced in the accretion disk (or the stream-impact region) 
we find that it is possible to create a consistent set of sys- 
tem parameters. As an X-ray burster, the compact object 
is a neutron star, whose mass we will assume is close to the 
canonical value of 1.4Mq. Furthermore, from our proposed 
short orbital period we can also infer that the donor star, if it 
is a main sequence star, must be of a late spectral type such 
as a M5V dwarf (which has a mass of 0.14Mq). Under these 
assumptions we can use Kepler's law to calculate that the 
binary radius is u—Q.QRq. Furthermore the donor's Roche 
lobe would then be 0.2_Rq (Paczynski 1971), i.e. comparable 
to the radius of such a late type star. 

Interestingly, de Jong et al. (1996) derived a mean disk 
opening angle for LMXBs of ~12°, making it possible that 
the disk in SerX-1 is indeed flared sufficiently to shield the 
donor star completely from irradiation. This could explain 
the absence of any indication of emission from the irradi- 
ated donor star in SerX-1. On the other hand, the presence 
of ionised Nitrogen, combined with the absence of Carbon 
suggest a more massive donor star. SerX-1 does not show 
any hint of the typically observed C HI A4647/4650 in its 
spectrum (see e.g. Steeghs & Casares 2002), while the curi- 
ous emission features in Fig.[T](those labelled with "?") could 
all be Nitrogen transitions. This strange abundance could be 
due to the CNO fusion reaction cycle (e.g. Shen & Bildsten 
2007), where proton capture on ^"^N is the slowest step. Since 
the CNO cycle is only dominant in stars more massive than 
our sun, this would imply that the donor in SerX-1 was orig- 
inally much more massive and that the processed material 
was mixed to the surface by convection. However, we know 
from Hynes et al. (2004) that Ha has a similar strength as is 
typically observed in LMXBs (note that from our Fig. [T] it is 
not even clear if H/3 is detected) , suggesting that the outer 
layers of the donor in SerX-1 have not been stripped. It is 
therefore unclear if the spectrum of SerX-1 can account for 
an evolved low mass donor. 

Despite the caveat mentioned above we can still con- 
sider whether the emission is coming from the full outer ac- 
cretion disk or a localized spot in the disk such as the stream- 
impact region using a low mass donor. In the first case the 
observed radial velocity curve should trace the motion of the 
compact object, while in the second case it would be that 
of the Keplerian orbit of the impact region. In both scenar- 
ios we need an estimate of the size of the accretion disk, 
and we assume that it reaches at least the truncation ra- 
dius (i?trunc). Using the system parameters described above 
we estimate that i?trunc=O.6Oa/(l-|-g)=O.457?0 (Frank et al. 
2002). Assuming a Keplerian accretion disk, the velocity at 
the edge of the disk would then be ~770 km s~^. If the 
line emission arises in the stream-impact region, which we 
expect to be close to the edge of the disk, an inclination 
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of i~3° is needed to explain the observed 38 km s~^ semi- 
amplitude. Furthermore, this inclination would also imply a 
non-projected FWHM of ~3650 km s~ . Since these num- 
bers, in particular the true FWHM, are rather extreme we 
think this unlikely. 

For the case of disk emission on the other hand, we 
can estimate an inclination of the system of i~9° using 
the mass function {f{M2)^M2sm^i/{l+qf=KlP/2-KG). 
This suggests that the emission lines would have a non- 
projected FWHM of ~1150 km s~^, which is compatible 
to the FWHMs observed in high-inclination LMXBs such as 
X1822-371 or GRMus (Casares et al. 2003; Barnes et al. 
2007). Furthermore, if SerX-1 is observed at such low in- 
clination the spectral resolution of our dataset would not 
allow us to resolve the double-peaked structure that is typ- 
ical for emission lines that are formed in the accretion disk 
(and have typical peak separations of ~600 km s~^), and 
the emission lines become indistinguishable from a single 
gaussian. 

Obviously, more complicated system geometries can be 
envisaged. One such complication, which is most likely im- 
portant in SerX-1, is the presence of strong tidal forces. For 
example Whitehurst (1988) showed that the accretion disk 
becomes tidally unstable and asymmetric for systems with 
an extreme mass-ratio (i.e. q^O.3). Since this is most likely 
the case in SerX-1 the accretion disk will not be a simple 
Keplerian one as we have assumed here. However, we do not 
think that this will affect our main conclusion, namely that 
we have discovered a ~2 hr periodic variability that we ten- 
tatively identify as the orbital period. If true, SerX-1 is a 
low inclination system with the spectral lines dominated by 
emission from the accretion disk. 
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